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ABSTRACT 

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array 
(ALMA) observations of the z=3.042 gravitationally lensed submillimeter galaxy HATLAS J090311.6+003906 
(SDP.81). These observations were carried out using a very extended configuration as part of Science Verifi¬ 
cation for the 2014 ALMA Long Baseline Campaign, with baselines of up to — 15 km. We present continuum 
imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), 
corresponding to an un-magnified spatial scale of — 180 pc at z=3.042. The ALMA images clearly show two 
main gravitational arc components of an Einstein ring, with emission tracing a radius of — 1.5". We also present 
imaging of CO 7=10-9, 7=8-7, 7=5-4 and H 2 O (2o2~ I 11 ) line emission. The CO emission, at an angular res¬ 
olution of-170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies 
between the CO transitions and compared to the dust continuum. Our detection of H 2 O line emission, using 
only the shortest baselines, provides the most resolved detection to date of thermal H 2 O emission in an extra- 
galactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure 
Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. 
Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the 
center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing 
galaxy. 

Subject headings: galaxies: high-redshift — submillimeter: galaxies — techniques: high angular resolu¬ 
tion — techniques: interferometric — galaxies: ISM — galaxies: individual: HATLAS 
J090311.6+003906 — gravitational lensing: strong 
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1. INTRODUCTION 

Dusty star-forming galaxies make up a significant fraction 
of the star formation rate (SFR) space density of the Uni¬ 
verse at z^2-3, where th ey are most numerous (e.g. Chapman 
120051: ICasev et sll 120141) . A powerful probe of their proper¬ 
ties is provided by gravitational lensing, with the magnifica¬ 
tion of both flux and apparent solid angle making it possi¬ 
ble to study in detail regions of luminosity-redshift space that 
would otherwise be inaccessible. Surveys that cover large ar¬ 
eas of sky, e.g. the Herschel Astrophysical Terahertz Large 
Area Survey (H-ATLAS ; lEales et al.ll2010|) . Herschel HER¬ 
MES (lOliver et al.l2012l) and the South Pole Telescope Survey 
d Vieira et al.l 11201 Of) have recently uncovered a large popula¬ 
tion o f strongly-lensed sources, c onfirming previous predic¬ 
tions (lBlainlll996t iNegrello et al.ll2007l). Other surveys, e.g. 
the Herschel Lensing Survey (lEgami et al.ll2QlQh are search- 
ing_for lensed sources in a targeted way. A review is given by 
iLutz et all (120141) . 

At submillimeter (submm) wavelengths, the study of red- 
shifted atomic and molecular lines in gravitationally lensed 
sources provides a useful diagnostic tracer of their interstellar 
medium (ISM) and star formation properties, active galactic 
nuclei (AGN), structure and dynamics (e.g. [S winban k et all 

2010. l20lH:lRiechers et al.ll201 lalilFu et aI1l2013l) . In addition 

to CO lines, the magnification provided by lensing also makes 
it possible to observe molecular species that would otherwise 

t, such as HCO + , HCN 

iRiechers et all 1201 lbl: 

hereafter OM13). H 2 0 
emission lines provide some of the strongest molecular 
lines in local Ultra-luminous Infrared Galax ies (U LIRGs) 
and starburst galaxies (e.g. iGonzalez-Alfonso et all (2004; 
Ivan der Werf et all 120101) . and strong water lines have now 
been detected in a number of high-z galaxies (OM13 and ref¬ 
erences therein). 

HATLAS J090311.6+003906 (hereafter SDP.81) is a gravi¬ 
tationally lensed submm galaxy (SMG) at z=3.042 that was 
detected in H-ATLAS by INegrello et al.1 ( 2010 . hereafter 
NE10). Based on spectroscopic redshifts from observations 
of several CO lines, NE10 determined that SDP.81 is a back¬ 
ground source at z=3.042 that is being lensed by a fore¬ 
groun d elliptical galaxy at z=0.299. Lens modeli ng analy¬ 
sis bv lDve et~aD (120141) and lBussmann et al.1 (I2013L hereafter 
BUI3) shows a magnification factor fi ~ 11. From NE14, the 
unmagnified infrared luminosity for the background source is 
5.1 x 1O 12 L 0 and the estimated intrinsic SFR is 527 M 0 yr _1 . 
SDP.81 has been studied in dust continuum emission using 
the Submillimeter Array (SMA) (e.g. NE10) and in molecu¬ 
lar lines (CO and H 2 0 (2 02 - In)) using the IRAM Plateau de 
Bure interferometer (PdBI; OM13), achieving angular reso¬ 
lutions of ~Q.6"-3". At this resolution, the dust continuum 
and molecular emission begin to resolve into multiple compo¬ 
nents, with the SMA and PdBI images indicating two gravita¬ 
tional arcs. 

Recently, the Atacama Large Millimeter/submillimeter Ar¬ 
ray (ALMA) carried out unprecedented high angular reso¬ 
lution continuum and spectral line observations of SDP.81 
as part of Science Verification (SV) of its 2014 Long Base¬ 
line Campaign. Details of the campaign, whose aim was to 
demonstrate the scientific capability of ALMA to observe on 
baselines of ~10 km or more, are described in an accompa¬ 
nying paper (IALMA Partnership et al.ll2015l. hereafter AL 15). 
In this paper, we present the multi-wavelength SDP. 81 ALMA 


SV observations. Section [2| describes the observations and 
data reduction. In Section [51 we present the data and initial 
results, and in Section 0 we discuss initial conclusions and 
the potential these rich data offer for detailed future analysis. 

Throughout this paper, we adopt a cosmology with 
Ho=61 kms -1 Mpc -1 , D m =0.32, Ua= 0.68 (iPlanck et al.1 
2014). 

2. OBSERVATIONS AND DATA REDUCTION 

Observations of SDR 81 in continuum, CO lines and an H 2 0 
line were carried out in 2014 October as part of the 2014 
ALMA Long Baseline Campaign using ALMA’s most ex¬ 
tended configuration to date. These data are publicly available 
from the ALMA Science PortaQ. The ALMA array consisted 
of 23 nominal long baseline antennas in a preliminary config¬ 
uration, plus a number of other antennas at shorter baselines 
(see AL15). During the periods of observation, between 22- 
36 antennas were in the array (with the number of available 
antennas varying with observing band) and baseline lengths 
ranged from ~15 m to 15 km. Only 10% of the baselines 
were shorter than 200 m. Observations were carried out in 
Band 4 (~2 mm), Band 6 (~1.3 mm) and Band 7 (~1.0 mm). 
Initial test observations in Band 3 (~3 mm) proved unfruitful, 
primarily due to the faintness of the continuum in this band 
due to the steep falloff of the thermal dust spectrum. 

For each of Bands 4, 6, and 7, the total available 7.5 GHz 
bandwidth was divided into four spectral windows (spws): 
three 15.6 MHz channel width spws for the 2.0, 1.3, and 
1.0 mm continuum, respectively, and one spw for the tar¬ 
geted spectral line. For Band 4, a 0.976 MHz channel width 
spw was centered on the redshifted CO 7=5-4 line (v res t = 
576.267 GHz); for Band 6, a 0.976 MHz or 1.95 MHz chan¬ 
nel width spw (online spectral averaging was applied for some 
executions) was centered on the redshifted low-excitation wa¬ 
ter line H 2 0 (2 02 - In) (v res t = 987.927 GHz (£ up = 101 K)); 
for Band 7, a 1.95 MHz channel width spw covered the red¬ 
shifted CO 7=10-9 line (v res t = 1151.985 GHz). In Band 6, 
one of the spectral windows included the CO 7=8-7 line (v res t 
= 921.799 GHz). Observed frequencies are given in TablesQ] 
and0 Other details of the observations, tunings and correla¬ 
tor setup are given in AL15 and Appendix lAl 

U v-coverage was critically important for imaging this tar¬ 
get, due to both its complex morphology and its location near 
the celestial equator that resulted in non-optimal horizontal 
uv- track£3- A relatively large amount of total observing time 
(^9-12 hours per band, or 9-12 executions of the scheduling 
block) was therefore necessary to achieve the good hour angle 
coverage needed for good image quality. The total on-source 
integration times were 5.9, 4.4 and 5.6 hours in Bands 4, 6, 
and 7, respectively. 

The data were reduced using the Common Astronomy Soft¬ 
ware Application package fCASAPI: iMcMullin et al.112007 ). 
The procedure is described in the scripts available on the 
ALMA Science Porta0. Imaging was carried out using a ro¬ 
bust = 1 weighting of the visibilities. Given the preliminary 
antenna configuration, there was a lack of wv-coverage for 
200-500 m baselines, and thus the use of robust weighting 
(robust= 1) was critical for achieving acceptable image quality 
(see Appendix [Aland AL15 for further details). None of the 

18 http://www.almascience.org 

19 http ://casaguides. nrao. edu/index .php ?title=ALMA2014_LB C_S VDATA 

20 http://casa.nrao.edu 

21 http://www.almascience.org 


be difficult to obser ve at high redshif 
or H 2 0 lines (e.j y Gao et al. 2007; 
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Fig. 1.— High resolution ALMA images of SDP.81: 2.0 mm (Band 4), 1.3 mm (Band 6) and 1.0 mm (Band 7) continuum. The smallest synthesized beam is 
31 x23 mas, for the Band 7 data. See Table[l]for other beam sizes. In addition to the gravitational arc structures, weak unresolved emission is also detected at a 
position that is coincident with the center of the foreground lensing galaxy. 


data were self-calibrated, since our attempts to perform self¬ 
calibration showed that there was insufficient signal-to-noise 
on the longest baseline antennas. The resulting synthesized 
beam sizes range from ^23-60 mas. 

For the spectral lines, it was necessary to taper the CO line 
uv data to 1000 kA in order to achieve good detections, re¬ 
sulting in coarser angular resolutions of ~170 mas. For the 
H 2 O data, it was necessary to taper the uv data further, to 200 
kA (^0.9" resolution). The CO data were binned spectrally 
into channels 21 kms -1 wide. For the narrower H 2 O line, we 
binned the data into 10.5 km s -1 channels. All the spectral line 
data were continuum subtracted and were imaged using rest 
frequencies corresponding to z=3.042. 

The resulting synthesized beam sizes and rms noise levels 
for the continuum and spectral line images are given in Ta¬ 
bles \l\ and [3j respectively. 

3. RESULTS 

3.1. Dust continuum emission 

In Figure |T] we present the high resolution continuum im¬ 
ages at 151, 236 and 290 GHz. The highest angular resolu¬ 
tion achieved, 23 mas, is an impressive increase of a factor 
of ^20-80 in angular resolution compared to previous SMA 
and PdBI observations of this source and corresponds to an 
unmagnified spatial scale of ~180 pc at z— 3.042 (or of the 
order of a few tens of parsecs in the source plane). Contin¬ 
uum images tapered to a similar resolution as the spectral line 
images (~170 mas) are presented in Figured 

Two main gravitational arc components of an Einstein ring 
are observed in all images: a larger eastern arc (E) that con¬ 
tains the highest peak intensity, and a smaller western arc 
(W). Several distinct surface brightness peaks are observed 
in each arc. The morphology is consistent with previous ob¬ 
servations at similar wavelengths with the PdBI and SMA 
(OM13, NE10). The locations of the arc components in the 
highest resolution (236 and 290 GHz) images lie on a ring 
of '■'4.55" radius that is consistent with values for the Ein¬ 
stein radius 0g= 1.52"-1.62" derived by NE10, BU13 and 
Dye et al. (2014) from SMA 880 p m or HST data. Mea¬ 
sured ALMA flux densities, both spatially integrated over all 
the arc components and for the major arc components (E and 
W) separately, are given in Table Q] These were computed 


within apertures that matched the spatial extent of the con¬ 
tinuum emission in each band using the tapered images. The 
tapered Band 6, Band 7 and combined Band 6&7 continuum 
images also show a suggestion of low S/N features that appear 
to trace out a more complete ring in addition to the two main 
arcs; due to the low significance of these features we did not 
include them in our flux measurement apertures. Peak fluxes 
were measured from the high resolution images. The ALMA 
2.0 mm (236 GHz) continuum flux density is in agreement 
with the PdBI 244 GHz value (2741.2 mJy; OM13) within 
the uncertainties. The 1.0, 1.14, 1.3 and 2.0 mm flux den¬ 
sities are consistent with the SED of the background source 
presented by NE14. 

3.1.1. Detection of continuum emission at the central position of 
the foreground lens 

In all three bands, we have detected continuum emis¬ 
sion located at the central position of the foreground lens 
(SDSS J090311.57+003906.5). This is visible in all the im¬ 
ages in Figure |T] and the 2 mm continuum tapered image 
in Figure [2 From the 236 GHz image, the centroid of the 
emission peak is at 09 /2 03 m llL57 +00°39 / 06".53 (J2000). 
The angular distance from this position to the eastern arc is 
1.5240.02", which is consistent with the Oe values deter¬ 
mined by previous authors. The emission is unresolved, even 
at 290 GHz, where the beam size is only 31 x23 mas. At the 
redshift of the lens (z=0.299; NE10), 23 mas corresponds to a 
spatial scale of ~ 100 pc. The detection of this central emis¬ 
sion suggests two main possibilities: a previously undetected 
AGN in the lens or an additional image of the background 
source that is predicted to arise close to the center of the lens¬ 
ing galaxy (e.g. lHezaveh et al.ll2015l and references therein). 
The ALMA continuum flux values at this positon differ by or¬ 
ders of magnitude from the SED predictions of both the lens 
and source in NE14, and we also note that we do not detect the 
central source in the CO or H 2 O line data. Although the ellip¬ 
tical galaxy is massive (NE10 find a mass within the Einstein 
radius of log(M^)=l 1.56 M 0 ), we find no evidence in the lit¬ 
erature of any AGN activity. However, the spectral index we 
derive from the ALMA 151, 236 and 290 GHz flux densities 
in Tableland the 1.4 GHz FIRST flux density (0.61 mJy) is 
-0.4940.13, which implies that the majority of the emission 
we have detected originates from the core of the foreground 
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Fig. 2.— ALMA images with wv-tapering to 1000 kA (CO lines and continuum) or 200 kA (H 2 O line). Top: CO J- 5-4, 8-7 and 10-9 velocity integrated intensity. 
Middle: 2.0 mm, 1.3 mm, and 1.0 mm continuum. Bottom: Band 6&7 spectral index, 1.14mm continuum (combined Band 6 & 7 data; see Appendix lAi. and 
H 2 O velocity integrated intensity. Beam sizes are ~170 mas, except for H 2 O which has a larger ~0.9 7/ beam (see Table[3]for spectral line details). 


elliptical galaxy. 

3.1.2. Spectral index 

The combined Band 6 and Band 7 spectral index image, ob¬ 
tained from data tapered to 1000 kA (see Section [3, is shown 
in Figured Pixels <4 a have been masked. The mean spectral 
index in the unmasked regions is 2.34±0.61, with the mean 
value measured in the W and E arcs comparable within the 
uncertainties (2.45±0.72 compared to 2.30±0.57). The range 
of values is consistent with dust spectral indices from 1.4-4. 
We note, however, that from the flux densities given in Ta¬ 
ble [U which represent the total emission region as opposed to 
the higher S/N regions defined by the spectral index image, 
the average spectral index is 1.8±0.2. 


3.2. CO J=10-9, J-8-7 and J=5-4 line emission 
3.2.1. Properties of the CO images 

CO 7=10-9, 7=8-7 and 7=5-4 velocity integrated intensity 
images are presented in Figure [2 The CO morphology is 
broadly consistent with the overall two-arc morphology seen 
in continuum, with a larger and higher surface brightness E 
arc and a smaller and generally lower surface brightness W 
arc. The symmetric distribution of the continuum peaks, how¬ 
ever, is not matched in CO, with the CO emission appearing 
more clumpy throughout the arc structures. The CO emission 
also appears to trace a somewhat more extended and less well- 
defined ring than the continuum at the same angular resolu¬ 
tion. For CO 7=5-4, the peak integrated intensities are com- 
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TABLE 1 

ALMA CONTINUUM PARAMETERS FOR SDP.81 


Band 

Frequency 

(GHz) 

Component 

S u pk 

[pjy beam -1 ] 

D" rms 

[p Jy] 

0b 

[mas (deg)] 

S„ 

[mJy] 

D" rms 

[p Jy] 

Otaper 

[mas (deg)] 

4 

151 

Total 


8 

56x50(18) 

3.5±0.1 

11 

134x118(61) 


(2.0 mm) 

E 

520± 8 



2.83=0.1 





W 

353±8 



0.73=0.1 



6 

236 

Total 


10 

39x30 (20) 

26.33=1.3 

21 

164x114(47) 


(1.3 mm) 

E 

1413=10 



19.33=1.3 





W 

1123=10 



7.0=b0.4 



7 

290 

Total 


9 

31x23(16) 

37.73=1.5 

21 

170x106 (43) 


(1.0 mm) 

E 

124d=9 



27.73=1.0 





W 

112±9 



10.0±0.2 



6&7 

262 

Total 




32.9±1.0 

15 

163x110(44) 


(1.14 mm) 

E 

8063=15 



23.5±0.5 





W 

627±15 



7.0 ±0.2 




Note. — Column (3) Spatial components: eastern arc (E), western arc (W) and the total over all the emission. (4) Peak flux density 
in the specified arc component. Values were measured in the high resolution images whose beam sizes, 0 b , are listed in column (6), 
except for the combined Band 6&7 data for which we used the wv-tapered image (see Section^. Uncertainties are lcr where a is the 
RMS noise level in the high resolution images, given in (5). For Band 6&7 the lcr uncertainty is the RMS noise value given i n (8) . 
(7) The integrated flux density over the E and W arcs separately, and the total emission over all the arc components (see Section [TH . 
Values are measured in the tapered images, whose beam sizes, 0 taper , are given in column (9). Uncertainties were calculated from 
the larger of either the rms variation of the flux density in apertures placed at source-free locations in the images or [no. independent 
beams] 0 ' 5 x3cr where a is the RMS noise level given in (8). Values do not include 5% absolute flux calibration uncertainty. (8) RMS 
noise level in the nv-tapered images. Synthesized beam sizes of the high resolution and tapered images, and the beam position angle, 
are given in (6) and (9), respectively. 


TABLE 2 

ALMA FLUX DENSITIES AT THE POSITION OF THE FOREGROUND LENS 


Band 

Frequency 

S v pk 


(GHz) 

{pSy beam -1 ) 

4 

151 (2.0 mm) 

60±8 

6 

236 (1.3 mm) 

50±10 

7 

290 (1.0 mm) 

43 ±9 


Note. — Frequency is the observed frequency. Peak fluxes were mea¬ 
sured in the high resolution images whose beam sizes are given in Table |T] 
Uncertainties are lcr where cr is the rms noise level given in Tableland do 
not include 5% absolute flux calibration uncertainty. 


parable in both arcs and the morphology bears more similarity 
to the continuum. The peak CO /= 10-9 and CO /=8-7 inten¬ 
sities, however, are respectively observed in the E and W arcs, 
and the CO /= 10-9 emission exhibits a different spatial mor¬ 
phology altogether, with several compact regions of emission 
along the arc structures. This may be due to differences in 
the spatial origin of the higher-/ CO emission, or simply due 
to the lower S/N in the CO J- 10-9 data. Notably, the spatial 
locations of peak intensity are not coincident between the CO 
transitions. 

Fitted spectral profiles, spatially integrated over the two arc 
components, are given in Figure [3] The spectra were pro¬ 
duced by spatially integrating the CO emission over a poly¬ 
gon region encompassing the E and W arcs as defined on 
CO integrated intensity images produced from all pixels with 
a S/N>3 in the channels containing emission. All parame¬ 
ters were derived from images that had first been smoothed 
to the same resolution (0.2"). The CO parameters derived 
from Gaussian fits are given in Table [51 For comparison, 
integrated intensity values for the E and W components are 
also given, and were measured directly on the integrated in¬ 
tensity images. For all measurements we do not include 
any low lev el em ission outside the main arc components 
(see Section 13.1b . The ALMA CO J- 5-4 integrated inten¬ 
sity measured over both the arc components, peak intensity 
and spectral shape are consistent with the PdBI CO J- 5-4 



-500 0 500 -500 0 500 

Velocity (km/s) Velocity (km/s) 


Fig. 3.— ALMA CO J= 5-4, CO 7=8-7, CO 7= 10-9 and H 2 0 (2 02 -ln) 
spectra for SDP.81. Spectra were spatially integrated over regions containing 
the E and W arc components (except for H 2 0, for which we only included E; 
see Section [Til in the nv-tapered images and the resulting spectra fitted with 
Gaussians. The two, redshifted and blue-shifted, components found for the 
CO transitions are marked with ‘R’ and ‘B’, respectively. 


results presented by OM13 (the PdBI integrated intensity is 
7.0±0.4 Jy kms -1 ). The three CO transitions have asymmet¬ 
rical line profiles (Figure[3]), with a red-shifted component (R) 
that is stronger than the blue-shifted component (B), which is 
again consistent with OM13. The B component is relatively 
brighter than R in the higher-/ CO lines, suggesting that the 
B emission originates from hotter gas than the R component. 
The line widths are in good a greement w ith the CO / = 1 -0 and 
CO J- 5-4 values found by Fraver et all (1201 lb and OM13, 
and arejvithin the range found for other typical SMGs (e.g. 
iGreve et aI1l2005 ). 
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TABLE 3 

ALMA CO and H 2 0 line parameters for SDP. 81 


Band 

Line 

^obs 

(GHz) 

Component 

Iipk 

(mJy beam -1 ) 

h 

(Jy kms -1 ) 

A Vi 

(km s -1 ) 

V 

(km s -1 ) 

nU 

(lo 8 L 0 ) 

/iL ; 7l0 lu 
(K km s -1 pc 2 ) 

6taper 

[mas (deg)] 

4 

CO 7=5-4 

142.570 

Total (fit) 


9.2±1.2 



9.6±0.3 

15.6±0.5 

155x121 (57) 




R 

26.5A0.2 

6.5±0.3 

231 ±6 

72±5 

6.8A0.3 

11.0A0.5 





B 

9.0A0.2 

2.7±0.3 

282±40 

-174±19 

2.8A0.3 

4.6A0.5 





E 


7.2±0.2 









W 


1.2±0.2 






6 

CO 7=8-7 

228.055 

Total (fit) 


9.3±0.4 



15.5A0.7 

6.2A0.3 

169 x 117(47) 




R 

28.0A0.2 

6.6±0.4 

221±11 

77±6 

11.0A0.7 

4.4A0.3 





B 

10.2A0.2 

2.7±0.4 

249±44 

-204±15 

4.5A0.7 

1.8A0.3 





E 


7.2±0.4 









W 


2.0±0.1 







H 2 0(2o2-ln) 

244.415 

Total (fit) 

9.4A0.9 

2.0±0.1 

196±15 

76±7 

3.6A0.2 

1.2A0.1 

852x721 (45) 




E 


1.6±0.1 









W 


0.15±0.06 






7 

CO 7=10-9 

285.004 

Total (fit) 


3.2±0.2 



6.7A0.1 

1.4A0.1 

172x112(42) 




R 

10.1 ±0.3 

1.8±0.2 

170±18 

75±8 

3.7A0.1 

0.8A0.1 





B 

6.7±0.3 

1.3±0.2 

188±30 

-204±12 

2.7A0.1 

0.6A0.1 





E 


1.2±0.1 









W 


0.4±0.1 







Note. — Parameters and errors for the CO and H 2 O lines are derived from Gaussian fits (Figure [3 for the two spectral components R and B (Section [±2l . 
All values are derived from the wv-tapered images, whose beam sizes and position angles are listed in the final column. v 0 b s is the observed frequency at z=3.042, 
which corresponds to zero velo city of fset. 7/ is the velocity-integrated intensity of the CO or H 2 O line and AV/ is the line width (FWHM). Luminosities, L/, are 
calculated following 1 Solomon et ali ( 1992 ). p is the magnification factor (Section |T}. Values of the int egrate d intensity for the two main spatial components (E 
and W arcs; Section [TT} are also given, measured directly on the integrated intensity images (Section [±2. IF Uncertainties on Iipk are lcr, where a is the rms 
noise level per channel. Rms noise levels for CO 7=5-4, CO 7=8-7 and CO 7= 10-9 are 0.20, 0.15 and 0.25 mJy per 21 km s -1 channel, respectively. For H 2 O, the 
rms noise level is 0.9 mJy per 10.5 km s -1 channel. Error values do not include 5% absolute flux calibration uncertainty. 


3.2.2. Velocity structure 

In Figure 0 we show the CO /=8-7 channel maps corre¬ 
sponding to the R component seen in Figure [3] in which a 
north-south velocity gradient is observed for the E arc. We 
do not show channel maps for the B component since orga¬ 
nized velocity structures are not clearly discernible in B. Sim¬ 
ilar properties are seen for the CO /=5-4 emission. This ve¬ 
locity structure of R is also demonstrated by CO /=8-7 inte¬ 
grated intensity and velocity field images produced separately 
for the R and B components (Figure [5]). Here, in addition to 
the N-S velocity gradient seen for R, it is also clear that the 
R and B image components have different spatial locations 
in the image plane. Spatially, the E arc corresponds predom¬ 
inantly to the R spectral component, with only a small part 
of the E arc that is blue-shifted. The emission in the W arc, 
conversely, has a symmetrical profile with roughly equal con¬ 
tributions from the B and R components. A velocity gradient 
is apparent for the W arc that encompasses both the B and 
R components, which can be seen in the velocity field image 
in Figure 0 The presence of two distinct spectral compo¬ 
nents, and the observed velocity structure in Figures 0] and 0 
is suggestive of complexity in the source structure, which has 
been previously suggested from lens modelling analysis (e.g. 
BU13). 

3.2.3. CO line ratios 

From the CO line luminosities given in Table 0 we 
find CO line brightness temperature ratios of rs5=0.5±0.1, 
rio8=0.2±0.1 and ri 0 5=0.1±0.L Taking the CO 7=1-0 value 
of 1.1 Jykms -1 from iFraver et al. ( 12011 ) we also find a ra¬ 
tio r5i=0.3= b0.1. This i s cons is tent with the valu e r3i=0.5 
reported bv IFraver et al.1 (120111) . iLupu et aD (1201 2b detected 
high-/ CO lines in SDR 81 but some were blended with other 
lines so we do not include them here. The values of the line 
ratios suggest the presence of a low-excitation gas compo¬ 
nent, which is consistent with results found for other SMGs 


(e.g. ICarilli et al.1 l2010h iDanielsonl 120 1 lb iHodge et aD 120 1 3b 
iBothwell et al.ll2013h and the suggestion that most SMGs con¬ 
tain a significant proportion of cool, moderate-density, ex¬ 
tended gas (e.g. llvison et ali 1 2011 ). Taking the total IR lu¬ 
minosity /jLir =5.4x 10 13 L 0 (uncorrected for magnification) 
from NE14, the ratio L IR fL' C0 is 346±62 L 0 (K kms -1 pc 2 ) -1 
for CO /=5-4, which istowards the upper end of the range for 
local ULIRGS (ISolomon et al.l fl997) and i n good agreement 
with the value found for typical SMGs (e.g. Greve et al. 2005, 
(360±140 L 0 (K km s -1 pc 2 ) -1 )). 

3.3. Water emission 

The H 2 O (2o2- 111) image is shown in Figure [2] Emission 
is clearly detected in the E arc and exhibits a morphology that 
appears consistent with those seen in continuum and CO. De¬ 
spite the tapering required to achieve a good detection, this 
is the highest resolution detection of thermal H 2 O emission 
in this source to date. For the W arc, there is a suggestion 
of possible emission, with integrated and peak intensities that 
are ~2cr in the integrated intensity image. The velocity range 
associated with this feature (^200 kms -1 ) is consistent with 
the B component. Due to the low significance of the W arc 
detection we do not include it in our spectrally fitted measure¬ 
ment. The ALMA integrated intensity and spectral profile are 
in good agreement with the PdBI value (1.8±0.5 Jy kms -1 ). 
The spectrum (Figure 0) shows that the H 2 O emission origi¬ 
nates from the R component. The central velocity of the H 2 O 
line profile agrees with all the CO line profiles within the un¬ 
certainties of the fits. 


3.3.1. H 2 O line ratios 

The H 2 O/high-/ CO line ratio is widely considered a 
useful diagnostic o_f_ PDR versus XDR conditions (e.g. 
[Gonzalez-Alfonso et al.11201 (1). We find an H 2 0/CO/=8-7 




































ALMA Long Baseline Observations of SDP.81 at z=3.042 


7 



Fig. 4.— Channel maps for the CO 7=8-7 emission, showing the red-shifted component, R. The R component was defined from -55 to 239 kms 1 based on 
the spectral profiles in Figure[3] Note that for clarity of the figure we do not plot the three outermost channel ranges. As a guide, contours of the combined Band 
6&7 continuum emission are shown in one panel, at 0.2 and 0.4 mJy beam -1 (15 and 30(j). Last panel: CO 7=8-7 velocity field, for both B and R components. 


line ratio, Ih 2 o/Ico of 0.2±0.1. This is relatively low com¬ 
pared to results in the literature for other active galax¬ 
ies. lOmont et~aD (120111) find a value of 0.5±0.2 for 
SDR 17b, while values of 0.4±0.2 and 0.75±0.23 have been 
found for the Cloverleaf and Mrk231 dBradford et all 120091: 
iGonz alez-Alf onso et al.ll2010h . Taking a value of L IR as in 
Section l±2.31 we find an L H 20 ^ir ratio of 6.7 x 10 -6 , which 
is similar to the value found for Arp 220 (OM13 and refer¬ 
ences therein). Both the intrinsic values of h H 2 o and h IR are 
at the lower end of the range for other high-z sources in the 
sample of OM13. 

4. DISCUSSION 

The results presented here on SDR 81 provide a first view 
of the wealth of information contained in these ALMA data. 
Providing a leap in angular resolution compared to previous 


observations, these data demonstrate the power of ALMA to 
image gravitationally lensed systems with high resolution and 
high fidelity. 

With the continuum images, we have presented strong de¬ 
tections of thermal dust emission from the z=3.042 back¬ 
ground source at unprecedented angular resolutions, with 
beam sizes as small as 23 mas corresponding to spatial scales 
in the source on the order of an impressive few tens of parsecs. 
At this high angular resolution we clearly detect the two main 
arc structures of the Einstein ring, but using only the shorter 
baselines we find evidence of low-level emission that traces a 
more complete ring. 

We have also detected weak (^5cr) continuum emission that 
is spatially coincident with the center of the foreground lens, 
and find a 1.4-290 GHz spectral index of ~-0.5 that implies 
we have detected emission from an AGN in the foreground 
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11.68 11.60 11.52 9:03:11.44 

Right Ascension (J2000) 


Fig. 5.— CO 7=8-7 integrated intensity (top) and velocity field (bottom) images for the B and R components. The B and R components were defined from 
-307 to -76 kms -1 and -55 to 239 kms -1 , respectively, based on their spectral profiles (Figure^- 


elliptical galaxy that must have a very low accretion rate. 

Our detection of H 2 0 (2 0 2 —In) in SDP.81, at 0.9" reso¬ 
lution, is the highest resolution detection to date of thermal 
water emission in an extragalactic source. H 2 0 emission from 
the W component is tentatively detected for the first time, with 
a contribution to the total H 2 0 emission of only a few per cent. 

We have also presented detections of three transitions of 
CO, J= 5-4, /=8-7, and J= 10-9, at 170 mas resolution, and 
shown that while they spectrally have similar red- and blue- 
shifted components, their spatial morphologies are rather dif¬ 
ferent. The properties of the CO and H 2 0 line data indicate 
that SDP.81 may be a complex source, _as was suggested by 
previous authors (e.g. BU13. lDve et al.lll2014l) . 

The overall similarity in the spatial and spectral distribu¬ 
tions of the CO and H 2 0 lines suggests that both species may 
have a common origin, although the sub-thermal (<Cl) CO 
line ratios and strong detections of high-/ CO lines could sug¬ 
gest the presence of an extended cool and moderate-density 
gas component (e.g. [Harris etaD 112010 ). Our CO line mea¬ 
surements suggest that the higher-/ CO lines are less lumi¬ 
nous, which may indicate that any AGN present would be 
subdominant, and thus that a complex mixture of PDRs and 
XDRs would be required to adequately describe the excitation 
in SDP.81. The H 2 0/C0(8-7) ratio is lower in SDP.81 than for 
other high-z sources in the literature (e.g. lOmont et al.ll201 ll 
OM13), which suggests that either we are missing significant 
extended H 2 0 emission or that the H 2 0 and CO excitation is 
relatively low. The former does not seem likely, given the sim¬ 
ilarity between the ALMA and PdBI results. It is possible that 


differential lensing (e.g. [Blainll 19991:ISerieanlll2012l) is playing 
a role in the relatively weak H 2 0 emission in the W compo¬ 
nent. Without detailed modeling, any effect from differential 
magnification cannot be quantified, and in a highly magnified 
system such as SDP.81, differential lensing could potentially 
lead to, for example, changes in the spectral indices, the CO 
ladder or line ratios. 

Future studies will be able to quantify this effect in SDR 81 
through detailed fo reground mass modelling ( Serjeant 2012; 
iBourne et al.l l2014). and the combination of this high angular 
resolution ALMA data with near-IR data will provide con¬ 
straints on the differential lensing between the stellar and dust 
emission (e.g. ICalanog et al.1120141) . Detailed lens modelling 
is beyond the scope of this work, but will be crucial for un¬ 
derstanding the full nature of these ALMA observations. 

This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2011.0.00016.SV. ALMA is a partner¬ 
ship of ESO (representing its member states), NSF (USA) 
and NINS (Japan), together with NRC (Canada), NSC and 
ASIA A (Taiwan), and KASI (Republic of Korea), in cooper¬ 
ation with the Republic of Chile. The Joint ALMA Obser¬ 
vatory is operated by ESO, AUI/NRAO and NAOJ. The Na¬ 
tional Radio Astronomy Observatory is a facility of the Na¬ 
tional Science Foundation operated under cooperative agree¬ 
ment by Associated Universities, Inc. This research made 
use of Astropy, a community-developed core Python package 
for Astronomy dAstropy Collaboration et alJl2013h . We thank 
Alastair Edge for helpful comments. 

Facilities: ALMA. 
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APPENDIX 

DATA REDUCTION DETAILS 

The phase center was 09^03 m lU.61 +00°39'06 // .7 (J2000). Precipitable water vapor (PWV) values at zenith varied in the 
ranges 0.6-3.2, 0.5-3.1 and 0.3-0.7 mm for Bands 4, 6, and 7, respectively, which is better than average weather conditions for the 
given Band. The Band 6 and Band 7 data were corrected with the most accurate antenna pad positions that were measured by the 
end of the campaign (see AL15). The phase calibrator, J0909+0121, was ~1.65° away from the target and was typically observed 
every 78 s for all three bands. Bandpass and flux calibration were performed using observations of the quasars J0825+0309 and 
J0854+2006 (or in a few cases J0750+1231 or J1058+0133), respectively, in each execution of the scheduling block (SB). For 
each dataset, between one and four antennas were flagged during the calibration process. 

Continuum images were integrated over a bandwidth of ^7 GHz and produced using multi-frequency synthesis, a robust-\ 
weighting of the visibilities, and nterms= 1. The use of robust weighting (robust= 1) was found to be critical for achieving 
acceptable image quality given the lack of wv-coverage for 200-500 m baselines due to the preliminary antenna configuration 
(since a limited number of antenna relocations were possible during the period of the campaign; AL15). The total number of 
antennas in the array depended on observing band, with fewest antennas in Band 4, due to fewer available Band 4 receivers 
on antennas on the shortest baselines. Multi-scale imaging (ICornwelll1120081 ) was used, with scales of 0, 5 and 15 for both the 
continuum and spectral line imaging (0, 5, 15 and 45 in Band 4). As the source lies sufficiently close to the phase center, the 
images have not been corrected for the effect of primary beam attenuation. This effect will be <5% in all bands, and thus smaller 
than the overall flux calibration uncertainty. 

We experimented with imaging the combined Band 6 and Band 7 uv data with CLEAN mode nterms- 2, which attempts to 
account for the spectral index of emission across the observed frequency band dRau & Corn welllHo 111) . The SNR of the full- 
resolution uv data was found to be too low for the algorithm to determine the spectral index, except around the strong peaks of 
emission, which resulted in a flux scale in the combined image that was not reliable. We were, however, able to make a combined 
(1.14 mm) image using Band 6&7 continuum images that were wv-tapered to match the resolution of the spectral line data, which 
yielded a robust spectral index image at ^170 mas resolution. 
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